Prediction of semi-metallic tetragonal Hf 2 03 and Zr 2 03 from first-principles 
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Abstract 

A tetragonal phase is predicted for Hf203 and Zr203 using density functional theory. Starting from 
atomic and unit cell relaxations of substoichiometric monoclinic HTO2 and ZrC>2, such tetragonal structures 
are only reached at zero temperature by introducing the oxygen vacancy pair with the lowest formation 
energy. The tetragonal Hf203 and Zr203 structures belong to space group P4m2 and are more stable than 
their corundum structure counterparts. These phases are semi-metallic, as confirmed through further Go Wo 
calculations. The carrier concentrations are estimated to be 1.77 x 10 21 cm -3 for both electrons and holes 
in tetragonal Hf 2 03, and 1.75 x 10 21 cm -3 for both electrons and holes in tetragonal Zr 2 03- The tetragonal 
Hf2C>3 phase is probably related to the low resistivity state of hafnia-based resistive random access memory 
(RRAM). 
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Hafnia (HfC^) and zirconia (ZrC^) are found in a number of important technological 
applications. 1, 2 In particular, hafnia has become a key component in sub-micrometer silicon MOS 
technology as the current choice of high permittivity dielectric layer. 3 In addition, it is also a 
promising candidate material for resistive random access memory (RRAM), which is one of the 
leading technologies for the next-generation non- volatile memory. 4-7 The core element of RRAM 
is a metal/insulator/metal capacitor which is subject to an electroforming process, where a high 
electric field (some MV/cm) is applied across the capacitor to create conduction paths in the insu- 
lating thin film, here named filaments. These filaments are of unknown composition or shape, and 
can be easily disturbed under electrical stress, leading to a memory effect. Knowing the composi- 
tion of the filaments is crucial to the understanding of RRAM's physics. Previous work on Ti0 2 
RRAM reveals that the conductive filament is possibly due to Ti n 2n -i Magneli phases, where the 
value of n is mostly 4 or 5. 8 For hafnia RRAM, the structure of the conductive filaments has not 
been reported, though it is widely accepted that the filaments are associated to an oxygen-deficient 
phased— Since the impact of electroforming is expected to occur in small and random patches of 
the capacitor, experimental investigation of the filaments suffers from great difficulty. 

In the present paper we employ first-principles density functional theory 11 (DFT) calculations 
to identify O-poor stable compositions of hafnium and zirconium oxides which may be reach- 
able from the room temperature normal pressure monoclinic Hf(Zr)02 [m-Hf(Zr)02] with the 
assistance of an external source of energy, possibly an applied electric field. The processing and 
operation of oxide-based RRAM (OXRRAM) stimulated our search of a conductive phase in these 
materials. However, our predictions are quite general and may have broader implications. 

For DFT calculations, plane-wave based Vienna Ab initio Simulation Package (VASP) 
program 12 was implemented, using projector-augmented-wave (PAW) pseudopotentials 13 with Hf 
hp, 5d and 6s (Zr 4s, 4p, Ad and 5s) electrons and O 2s and 2p electrons in the valence. Generalized 
gradient approximation (GGA) was used for the exchange-correlation energy, within the Perdew- 
Burke-Ernzerhof (PBE) functional. 14 The plane wave energy cutoff is chosen as 500 eV which 
converges for all the involved compounds, and sufficiently dense Monkhorst-Pack 21 or T-centered 
fc-mesh was utilized for sampling the Brillouin zone. 

Because DFT/GGA usually underestimates band gaps, and may even deem a material metal- 
lic rather than semiconductor as the case of bulk germanium, 15 for the metallic or semi-metallic 
candidates we calculated the first order energy shifts with GW approximation 16 (GoWo) us- 
ing the ABINIT— program. Hf (Zr) semi-core electrons were explicitly included through 
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the 5s 2 5p 6 5d 2 6s 2 (4s 2 Ap 6 Ad 2 '5s 2 ) configuration while core electrons were replaced by Troullier- 
Martins pseudopotentials. 19 Convergence was achieved with 360 bands and a 15 Ha cutoff for 
the wave functions employed in the evaluation of the dielectric function and the E function. A 
plasmon-pole approximation 2 - was used. 

The m-Hf(Zr)0 2 unit cells were fully relaxed until all Hellmann-Feynman forces were less 
than 0.01 eV/A and all stresses were less than 400 MPa. The relaxed m-Hf(Zr)0 2 unit cell pa- 
rameters are a=5. 146 (5.219) A, b/a=l. 010 (1.012), c/a= 1.036 (1.036) and = 99.68° (99.68°), 
close to experimental values:^ a=5. 117(5. 151) A, b/a=1.011 (1.010), c/a= 1.034 (1.032) and /3 = 
99.22° (99.20°). The formation enthalpy of m-Hf(Zr)0 2 calculated with respect to h.c.p. Hf (Zr) 
and an isolated 2 molecule is -1166(-1106)kJ/mol, after adopting the 1.36eV energy correction 
for the 2 molecule provided by Wang et alP These results are in accord with the experimental 
values, which are -1 145 kJ/mol for hafnia and -1101 kJ/mol for zirconia. 23 

For reference we first calculated the formation energies of single and double oxygen vacancies 
in a 96-atom m-Hf(Zr)0 2 supercell (2 x 2 x 2). Aiming at energetically favorable metallic phases, 
only the neutral oxygen vacancy was considered. There are two inequivalent O sites regarding O 
coordination, namely the 3- and 4-coordinated O(A) and O(B) sites, respectively. The formation 
energy of a neutral oxygen vacancy is defined as 

-Efbr = Ed — Eq + /i (1) 

where is the energy of the defective supercell, E is the energy of the defect-free supercell 
and no is the chemical potential of oxygen. Under oxygen-rich condition /^o is usually set as 
one half of that of an 2 molecule. In this case the calculated Ef OT for a neutral O(A) vacancy in 

All inequivalent 
di-oxygen-vacancy sites 

(A1,A2) (A1,B3) 
(A1,A3) (A2,B2) 
(A1,A4) (Bl, B2) 
(Al, Bl) (Bl, B3) 
(Al, B2) (B2, B3) 

FIG. 1 : Unit cell used for calculations of the energy of formation of oxygen vacancies in Me0 2 (Me = Hf 
or Zr) with all metal and oxygen sites identified. The ten inequivalent pairs of oxygen sites are also listed. 
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FIG. 2: Density of states for (a) m-Hf(Zr)0 2 ; (b) Hf(Zr) 4 7 with one 0(A) vacancy; (d) Hf(Zr) 4 7 with 
one 0(B) vacancy. Black curves are for Zr-based compounds while red curves are for Hf-based compounds. 
All figures are aligned horizontally with respect to the O 2s band (not shown). The highest occupied 
molecular orbital levels are indicated by vertical dashed lines. 

m-Hf(Zr)02 is 7.13 (6.62) eV, while for a neutral 0(B) vacancy it is 7.00 (6.53) eV. The difference 
between the two formation energies, 0.13 (0.09) eV, is similar to Zheng et ah, 26 but larger than 
Foster et al. who reported a 0.02 eV difference in both cases.—'— Next we calculated the formation 
energies of di-oxygen-vacancy pairs in m-Hf(Zr)02. To this end, the 8 oxygen sites in a unit cell 
were named A1-A4 and B 1-B4, as shown in Fig. 1 . The distance between two di-oxygen-vacancy 
pairs is around 10 A, casting them as isolated pairs. Since m-Hf(Zr)02 possesses the baddeleyite 
structure with space group P2i/c, there are 10 inequivalent di-oxygen-vacancy pairs. The most 
energetically favorable pair is the (Bl, B2) pair, whose formation energy per vacancy is the same 
as of a single O(B) vacancy. 

Under usual experimental conditions the dielectric is placed between two metal electrodes 
where oxygen can migrate as an interstitial. We thus calculated the incorporation energies of 
oxygen, starting from its molecular form, into bulk h.c.p. Ti and f.c.c. Pt, two commonly used 
electrodes. The results are -6.24 eV and 0.91 eV, respectively. The formation energies of single 
and double oxygen vacancy in all three environments are thus compared in Tab. I. Notice that the 
formation energy of a neutral 0(B) vacancy in m-Hf(Zr)0 2 plus an oxygen interstitial in h.c.p. Ti 
is merely 0.76 (0.29) eV, which is attributed to the strong Ti-0 bonding. 
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The search for metallic states began with the Hf(Zr) 4 07 models generated by introducing one 
oxygen vacancy per 12-atom m-Hf(Zr)02 for two inequivalent cases, i.e., 0(A) and 0(B). The 
defective unit cells were fully relaxed until all Hellmann-Feynman forces were less than 0.01 
eV/A. Atomic coordinates, cell dimensions and shape were subject to relaxation. The formation 
energies per vacancy of the two Hf(Zr) 4 7 phases are almost the same as in the single oxygen 
vacancy cases. The resulting DOS are shown in Figs. 2(b) and 2(c). In either case, a fully occupied 
defect-induced band emerges in the band gap, indicating a semiconductor. While the DOS of this 
band is broad and high, it cannot account for the measured metallic state in hafnia-based RRAM, 
because in that state the resistance is on the order of hundreds of Ohms for 10 nm thin films. 29 

Chemical Vacancy Supercell Formation energy (eV/vacancy) 
formula site(s) units with oxygen going to 
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TABLE I: Formation energies of isolated single oxygen vacancies and di-oxygen-vacancy pairs, of isolated 
oxygen vacancy chains, and of substoichiometric Hf(Zr) 4 7 and Hf(Zr)20 3 in different structures derived 
from monoclinic Hf0 2 and Zr0 2 . In the table, Me stands for Hf or Zr; the parenthesed formation energy 
values are Zr-based while others are Hf-based. 
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However, the trend does hint that stronger off- stoichiometric hafnia or zirconia might undergo a 
phase transition from dielectric to metal. 

Hence, several Hf(Zr) 2 3 models were set up with two oxygen vacancies per 12-atom m- 
Hf(Zr)0 2 unit cell. For all inequivalent cases (Fig. 1) the cells remain monoclinic during re- 
laxation of the unit cell vectors, except for the (Bl, B2) case which suffers from a monoclinic-to- 
tetragonal transition. The tetragonal Hf(Zr) 2 03 [t-Hf(Zr) 2 3 ] phase (Fig. 3; structural parameters 
in Tab. II) is the ground state of all ten Hf(Zr) 2 3 candidates. It belongs to the D 2 d point group 
and PAm2 (No. 115) space group 31 . The Hf(Zr) coordination number is 7 as in m-Hf(Zr)0 2 , 
while 2/3 of the oxygen sites have coordination number 5 [named 0(A)] and 1/3 have coordi- 
nation number 4 [named 0(B)]. Symmetry analysis indicates that Hf(Zr) and 0(A) are at the 2g 
position while 0(B) is at the lc position 32 . The average Hf(Zr)-0(A) and Hf(Zr)-0(B) bond 
lengths are 2.295 (2.322) A and 2.089 (2. 134) A, respectively, compared with 2.084 (2. 1 17) A and 
2.209 (2.240) A in m-Hf(Zr)0 2 . Bader analysis reveals that the Hf (Zr) charge changes from 
2.73e(2.57e) in m-Hf(Zr)0 2 to 2.10e(2.02e) in t-Hf(Zr) 2 3 ; the 0(A) charge changes from - 
1.34e(-1.25e) to -1.39e(-1.34e); and the 0(B) charge changes from -1.39e(-1.31e) to -1.41e(- 
1.37e). 

Figure 4 shows the band diagram and orbital-projected DOS of t-Hf 2 3 . Results for i-Zr 2 3 
(not shown) are similar. The high symmetry points in the Brillouin zone are named according to the 
Bilbao crystallographic server. 31 A semi-metallic behavior is revealed in the overlap of the partially 
occupied valence band top and conduction band bottom located at different high symmetry points, 
R and Z. Both band edges are mostly derived from Hf 5d states. The semi-metallic character 
of the two compounds was confirmed by GoWo calculations of the many-body correction to the 
DFT/GGA energy levels. Figure 4 implies that the energy shifts are of the order of a few tenths 
of eV, retaining the semi-metallic nature of t-Hf 2 3 . The energy shifts obtained for Zr 2 3 (not 
shown) are similar. From the calculated band structures the densities of conduction electrons 
and holes can be obtained by integrating the electron occupation of the blue and orange bands in 

a c Hf z 0(A) Z Bulk modulus 

Hf 2 3 3.135A 5.646A 0.2553 0.1351 246GPa 
Zr 2 3 3. 174 A 5.763 A 0.2525 0.1367 228 GPa 

TABLE II: Calculated structural parameters of tetragonal Hf 2 3 and Zr 2 3 . 
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FIG. 3: Tetragonal Hf2C>3 and Zr203.— (a) Primitive cell with 5 atoms; (b) view of the structure along a- 
axis (upper) and c-axis (lower); (c) simulated powder X-ray diffraction patterns of Hf2C>3; (d) Fermi surface 
of Hf2C>3 at T=0 K. Big green balls: Hf or Zr; small red balls: O. 




-2 J 1 1 1 1 1 1 1 

X T Z R AM Oxygen Hafnium 



FIG. 4: Tetragonal Ett^Os electronic band structure (left) and orbital projected DOS (right). DFT/GGA 
(solid lines) and Go Wo (discrete marks) band structures are superimposed and aligned by their Fermi levels. 
Results for Zr20^ (not shown) are similar. 
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Fig. 4, respectively. The electron and hole densities are both 1.77 x 10 21 cm~ 3 for t-Hf 2 3 , and 
both 1.75 x 10 21 cm -3 for t-Zr 2 3 , typical of semi-metals. 

To evaluate the relative stability of this structure, we calculated the molar formation enthalpy 
of various Hf 2 3 , Zr 2 3 and Ti 2 3 models with respect to their corresponding metals and 2 . 
Still, the 1.36 eV energy correction to 2 was applied to all cases. The formation enthalpy of 
t-Hf(Zr) 2 3 is - 1700 (-1 666) kJ/mol, more favorable than fully relaxed corundum Hf(Zr) 2 3 , - 
1586 (-1580) kJ/mol. Nevertheless, a fully relaxed Ti 2 3 arranged in the PAm2 tetragonal struc- 
ture, possesses a formation enthalpy of -1576kJ/mol, less favorable than corundum Ti 2 3 whose 
formation enthalpy is -1598kJ/mol. These data confirm that for Hf(Zr) 2 3 the tetragonal P\m2 
structure is preferred, while for Ti 2 3 the corundum structure is preferred. 

To our best knowledge, the proposed t-Hf 2 3 and t-Zr 2 3 structures have not been reported 
before, though some published data may suggest their existance. Hildebrandt et al. 33 performed 
high-resolution transmission electron microscopy of a conducting Hf0 2 _ x thin film where the 
enlarged inverse Fourier-transformed images show a similar structure as in Fig. 3(b). Manory et 
ah— discovered two unidentified X-ray diffraction (XRD) peaks at 29 = 40° and 29 = 52° in 
hafnia films grown by ion beam assisted deposition at a transport ratio of 5 and an ion energy 
of 20keV. They attributed these peaks to a new tetragonal structure and suggested the Hf 2 3 
stoichiometry. However, they simulated their data with a PA/mmm phase with lattice parameters 
a=5.055 A and c=5.1 1 1 A, resulting in two small peaks around 40°. We calculated powder XRD 
patterns 35 for t-Hf 2 3 [Fig. 3(d)] using Cu Ka radiation (A = 1.5418 A), and found a (110) peak 
at 40.7° and a (1 12) peak at 52.5°, similar to data. Similar calculation for t-Zr 2 3 yielded a (1 10) 
peak at 40.2° and a (112) peak at 51.6°. 

In conclusion, we have predicted tetragonal semi-metallic Hf 2 3 and Zr 2 3 structures as the 
ground state highly oxygen deficient hafnia and zirconia which undergo a monoclinic-to-tetragonal 
phase transition. Their semi-metallic properties are characterized by an overlap of the valence 
band maximum and conduction band minimum at different points of the Brillouin zone, and by 
the low density of conduction electrons and holes. Also, i-Hf 2 3 may be the physical origin of 
the conductive state in hafnium-based RRAM. 
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